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Abstract 

We develop the extended dynamical mean field theory (E-DMFT) with a 
view towards realistic applications. 1) We introduce an intuitive derivation 
of the E-DMFT formalism. By identifying the Hartree contributions before 
the E-DMFT treatment, it allows to handle systems in symmetry breaking 
phases within a simple formalism. 2) We make a new implementation of 
E-DMFT through real Hubbard-Stratonovich transformation to decouple the 
non-local two-particle interactions. We apply it to a 3D U-V model, with U 
the on-site and V the nearest neighbor interactions, and investigate the be- 
havior of the various Green's functions, especially the density susceptibility, 
as the density instability is approached. We obtain the phase diagram at a 
finite temperature. 3) We present a formalism incorporating E-DMFT with 
Cellular DMFT. 4) We suggest an improvement of the E-DMFT approach 
by combining it with a generalized GW method. The method combines the 
local self-energy from E-DMFT and the non-local ones from the perturbative 
calculation of GW. We apply the method to a ID U-V model with two sub- 
lattices carrying different chemical potentials. By comparing with those from 
Density Matrix Renormalization Group (DMRG) calculations, we show the 
results are shifted in the correct direction due to the GW contributions. 5) In 
order to handle the generic Coulomb repulsion within E-DMFT, we describe 



a method to tailor E-DMFT so that proper momentum dependence can be 

kept in general response functions. 
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I. INTRODUCTION 



The Dynamical Mean Field Theory (DMFT) has been a powerful tool for the study 
of strongly correlated electronic systems. It has allowed us to gain new insights into non- 
perturbative problems such as the Mott paramagnetic metal to paramagnetic insulator tran- 
sition at finite temperatureJil. There are many current attempts to extend the scope of 
the DMFT approach in several directions: to include momentum dependence in the self- 
energyEHi, to account for the effects of finite range interaction^Hli, and to combine DMFT 
with realistic band structureEHli. 



The current paper is devoted to the investigation of the Extended Dynamical Mean Field 
Theory (E-DMFT), an extension of the original DMFT in order to take into account the spa- 
tially non-local interactions beyond the Hartree level. The idea of E-DMFT was developed 
independently in the studies of spin glasslill, systems with non-local Coulomb interactionillll, 
and the heavy fermion syste m0. The derivation of E-DMFT based on Baym-Kadanoff func- 
tional has been achieved in Ref.H. 



We present here several methodological developments which build on the E-DMFT ap- 
proach with a view to obtain a more realistic description of solids. Our goals are to describe 
a) the frequency dependence of the effective interaction and its effects on the single particle 
spectra, b) the effects of short range correlations, and c) a momentum dependent polar- 
ization. These three effects certainly present in realistic models of solidsM. We discuss 
them in this paper in the framework of model Hamiltonians in sections |VI1] , and |VII1 . 
It should be pointed out that the problem is not touched in this paper as how to obtain the 
parameters of the model Hamiltonian from the first principle calculations. 



In addition we present several technical advances for the analysis of E-DMFT equations. 
In sections |I| and |T|, we present a simplified derivation of these equations in a broken 
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symmetry phase. A method for handhng arbitrary interactions within E-DMFT-QMC by 
an interaction shiftH is also discussed in the section. We show the formahsm combining 
E-DMFT with cellular DMFT (C-DMFT) in §|V| and with GW approximationB in §0. 

To illustrate the ideas and the working of the methods we apply them to two problems. 
The transitions between a Mott insulator (MI), a band insulator (BI), and a Fermi liquid 
(FL) in a 3D U-V model is discussed in The model describes an electron system with 
an on-site repulsion U and a nearest neighbor interaction V. It was treated in simple DMFT 
at hall and quarte ri fillings and is relevant to the materials with charge ordered phasecBca. 
The problem of the transition between MI and BI phases in a ID U-V model with alternating 
chemical potential was discussed in the context of mixed stack organic compoundSS and 
ferroelectric perovskite^. We exhibit in §[VII| the implementation of E-DMFT-GW method 



on the ID model. 5pC| is the conclusion. 



II. MODEL HAMILTONIAN 

We start with the following Hamiltonian: 

1 3 

i,j o" i a=0 

1 3 ^ 

+ UY^ hi^hii + SiaViajpSjp. (1) 

i i,j a,/3=0 

The operators Sia *= cl„T^^,Cia' with t°' being the Pauli matrices for a = 1,2,3 and 
■7"° = -^2x2, the identity. So the zero-th component of Sia is the charge density and the 
rest the spin operators. Similarly, hia=o represents the chemical potential and the other 
three, with a = 1, 2, 3, the external magnetic field. 
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In Eq.(|T]) the on-site part of the interaction is isolated in U and the off-site parts are 
described by ^10,^73, hence Via^ip = 0. We have allowed all the possible forms of the instan- 
taneous direct and exchange interactions, but excluded some others, like the pair hoppings. 
Unless otherwise specified, the form of the interaction Viajf3 will be generic with only the 
simple requirements that the interaction be translational invariant, Vi^jf^ = Va^pdi — j\), 
and symmetric. Via j 13 = VjPM- -Due to the possible ionic screening and the super-exchange 
mechanisms, we have the freedom to certain extent in choosing U and the components of V 
independently at the level of the model Hamiltonian. While it is reasonable to take f/ > 0, 
in the off-site interaction matrix V all possibilities are allowed if some effective forms of 
the interactions are under consideration. In other words, the matrix V is not necessarily 
positive- or negative- definite. 

In this paper we keep for simplicity h^^h^ 7^ and set hi,h2 = 0. This choice is more 
general than it appears since we can always choose the direction of the effective magnetic 
field to be the z-direction by appropriately rotating the system. Furthermore, we assume 
that the spin symmetry breaking, if it may happen, is also along the z-direction. This 
assumption is very natural if Viaj/3 is symmetric under the permutation among the three 
directions. However, it is actually more restrictive since it requires that if Va,i3{\i — j\) 
produces any easy-axis, the direction of the easy-axis should be also in the z-direction. We 
can then define a spin-dependent chemical potential: 

Before deriving the E-DMFT formalism, we need first to separate out the Hartree con- 
tributions from the interactions. It is known from the Baym-Kadanoff functional analysis 
of E-DMFT0 that in the phases with broken symmetry, Hartree terms contribute to the E- 
DMFT by shifting the chemical potential^. In the current case, due to the spin exchanges, 
this shift is spin-dependent. It is also known that in the circumstance of C-DMFTii, one 
needs to handle separately the Hartree contributions from the non-local interaction across 
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the cluster boundary. With all these motivations, we rewrite the Hamiltonian as follows: 

H = EE(^ij4ci<x + h.C.) -Y.^€n^a 

i,j f i,a 

1 3 ^ 

+ UY, ■ ■■ nil : +-Y^ E : Sia ■ Vi^jfi : Sjp : (2) 

i i^j a,f3=0 

where the normal ordering of the operators is defined by : O : =^ O — {O) with the average 
over the ground state. The effective chemical potential is defined to be 

We work with the functional integral representation of the partition function at finite tem- 
perature < T = We want to perform a Hubbard-Stratonovich transform and 
decouple the \^-interactionsi. Since ultimately we will map the many-body system to an 
impurity model and solve the self-consistent impurity problem via numerical techniques like 
QMC, it is desired that the Hubbard-Stratonovich transform be reaS To this end, we need 
first to add the identity matrices Xaphj to the off-site interaction to ensure that, for given 
a and /?, the matrix 

is positive-definite. The minus sign in front of the bare interaction in the above equation 
is needed for a real Hubbard-Stratonovich transform. Practically we can take any value 
of Aq/3 as long as it is greater than the biggest eigen value of the matrix \Vap\ whose el- 
ements are defined by [Kt/?]*^ *= Viajp- To keep the symmetry, we require Xap = X/sa- 
Using Hubbard-Stratonovich transform, one can write the partition function in terms of the 
following functional integral: 

Z = J P[4(r), Q,<,(r); 0,,,(r)] exp(-5) (5) 

with the Euclidean action: 
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Jo T:" r1 ^ 



+ U^^j:--Mr) ■■■■ n^ir) ■■ E <PUr)[Vap]^%Ar) ^T^T. M^) S^a{r) :}, (6) 

i i,j a, (3=0 i a=0 

where 

= U + Vio.iO ~ ~ ^2,i2 ~ ^i3,j3 = U + Aqo — An — A22 — A33 

While the A dependences are exphcitly contained in the effective on-site interaction U'^^ 
and the interaction matrices Vap-, they cancel each other exactly so that all the physically 
measurable quantities are independent of A. The electron-phonon vertex defined in Eq.(^ 
is local. In an insulator with a long range Coulomb interaction, one needs a different choice 



of the vertex, as will be discussed later in 5 VIII 



Several remarks are in place. First, since our strategy is to separate out all the Hartree 
contributions from the interactions, we certainly need to treat the auxiliary field in the same 
way. However, due to the identity 

{^^o) = TY.^l[yo.pU■■S,p :)=0, 

j 13=0 

we need not bother to normal order the fields here. The current formalism of E-DMFT 
can be extended directly to a general bose fieldi which can develop non-zero expectation 
values, if we replace 0jQ, by : (pi^ :. Second, there is an arbitrariness in choosing the sign of 
the electron-phonon vertex which is reflected in the "±" signs appeared in Eq.(|^). It comes 
from the freedom in the Ising {Z2) symmetry of the order parameters. For definiteness, we 
take '— ' from now on. Finally, one can also try to decouple the on-site interaction by using 
yet another continuous auxiliary field. Because this interaction does not get any renormal- 
ization as we derive the E-DMFT, and because for the local density-density interaction there 
is the more efficient Hirsch-Fye algorithmil which uses the discrete Ising auxiliary field, we 
will not implement this here^. 
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For future reference, we first define the following Matsubara Green's function: 

G..,v(r - t') _(T,c,.(r)4,(r')), (7) 

Xiajfiir - r') = -(T^ : Si^ir) :: Sjpir') :), (8) 

Diajpir - t') -{Tr : (^Ur) :: 0,^(^0 ■)■ (9) 

Notice that the auxiliary field Green's function D depends on the value of A introduced to 
make the effective interaction matrix V positive definite while the Green's functions for the 
electrons do not. The auxiliary phonon Green's functions are related to the electronic two- 
particle response functions by the following identity, which can be derived by integrating 
out the auxiliary degrees of freedom in the phonon Green's function, 

X-\k,lUJn) = Vk + Il-'{k,lUJn) (10) 

with the self-energy given by the Dyson equation for the phonons: 

U{k,lUJn) = -V,-' - D-\k,lUJn), (11) 

where cUn = (27r//5)n for integer n. The above two equations are in the matrix form in the 
charge-spin space labeled by a = 0, 1, 2, 3. The bare interaction vertex —V^ plays the role as 
the free propagator and is defined as [Vklap i^/^^) J2j,iVja,w • {j — I)]. For later 

use, we also write down the Dyson equation for the electron Green's function: 

E{k,ipn) = ipn-tk- G~\k,ipn). (12) 

with Pn = (27r//9)(n -|- 1/2) for integer n. The equation is in the 2x2 matrix form labeled 
by the spin a J,. 



With the above preparation, we now derive E-DMFT. We will consider two situations. 
In the next section we describe the first case, the E-DMFT for a homogeneous system in 
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which all the lattice sites are equivalent to each other. In this case E-DMFT approxima- 
tion amounts to integrating out all but one lattice site to get an effective single site action 
which is equivalent to a self-consistent impurity model. All the sites that are integrated out 
contribute to self-consistent baths of free fermions and bosons for the impurity model. Here 
we should have a homogeneous chemical potential = fJ^a- In the second case (see §|Vp we 
explore the E-DMFT for systems with two mutually penetrating sublattices. We present a 
formalism combining E-DMFT and C-DMFT with a cluster of two neighboring sites. While 
it is natural to apply this formalism to a system with two non-equivalent sublattices, it can 
also be applied to the homogeneous situation. In the latter case, the purpose is to treat 
more accurately the spatial correlations than the single site E-DMFT. This combination of 
E-DMFT and C-DMFT can be formulated for clusters of arbitrary size. 



III. E-DMFT OF A HOMOGENEOUS SYSTEM 

We first explore the single site E-DMFT applicable to a homogeneous system. By keeping 
one lattice site while integrating out the rest using the cavity constructionS, we obtain the 
E-DMFT effective actioni 

Sf = - fdr f dr'Y.clAr)g-\r-T')c,Ar') 
Jo Jo „ 



-\ f dr f dr' Y: ■■ Mr) : p-j(r - r') : ^r') 

^ Jo Jo „^n 



+ [/^fif fdr : no,T(r) :: no,i{r) -.- f drj^: 0oa(r) :: So^r) ■ (13) 

-'0 Jo t^r, 



dr 

0=0 

This is an impurity model with both electron and (auxiliary) phonon degrees of freedom. 
The impurity model Green's functions are identified with the local Green's functions of the 
lattice model, namely [see Eqs.(|^)-(^)], 
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GlJ^Mr - r') = Go.,o.'(r - r'), " = Xoa.o^r - r'), I^L°/3'(r - r') = DoaM^ " r'). 

(14) 

The E-DMFT self-consistent loop is formed as follows. One starts with the effective action 
Eg . (p!3D and measures the local electron and phonon Green's functions as define through 
Eq.(|l4D. Then one calculates the self-energies using the local version of the Dyson equations: 

^AWn) = Q-\iPn) - [G':r\Wn), (15) 

n(ia;„) = V-\iuJr,) - (16) 



One of the basic assumptions of DMFT is the locality of the self-energyu. Notice that 
from Eq. (|To|) the phonon self-energy is different from, although related to, the response 
functions and the locality of the former does not imply the locality of the latter. Under this 
assumption, the same local quantities can be calculated, by using Eqs. (|T0|)-([T2l), 

G':\^Vn) = E G.{k, ipn) = -ik- fif - ^A^Pn)]-\ (17) 

k k 

X'°'{^u;n) =Y.x{k,iuj^) = Y.[Vk + Ii-\iuj^)]-\ (18) 

k k 

D'°\lUn) =J2D{k,lUn) = -Ei^fc"' +n(^Wn)]-' (19) 
k k 

In the DMFT with phononS, one substitutes the results of Eqs.(p!7D and (|19D back in the 
local Dyson equations, (|I5D and (0), and generates a new set of the Dynamical Weiss func- 
tions Q and V which serve as the starting point of the next iteration. This completes the 
self- consistent loop. 

It is interesting to notice, however, that one has another choice to form the self-consistent 
loop through the electronic two particle Green's function instead of D^"'^, as used in 
Refs. PIlslO] . It is important that the two different procedures are compatible with each 
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other. By combining Eqs. ([ISD and (|T9D , one can find the following relation in the matrix 
form: 



[x'^T'i^^n) = -V{lUJn) + Tr\lUn) (20) 

where the local phonon Dyson equation has been used. This is nothing but the iden- 
tity one can derive directly from the effective action Eq.([l3D and is the local version of the 
electron-phonon identity, Eq. (|lO|) . Hence the two routines, which correspond to making the 
DMFT approximation at the different stages of the formulation, one directly from the non- 
local electron interaction and the other after the Hubbard-Stratonovich decomposition, are 
indeed equivalent, as long as the same interaction vertex defined in Eq.(|) is used. 

There is, though, still one point left out, that is the arbitrariness of the value of the con- 
stants Xap in Eq.(^. They were introduced in our formalism for the purpose of ensuring the 
applicability of the real Hubbard-Stratonovich transform and are thus unphysical. While it 
is obvious that the artificial effects from the A dependences cancel out exactly in the original 
action, Eq.(|^), it is not clear if the same thing happens after the E-DMFT approximation. 
On the other hand, the earlier versions of E-DMFTi0i use only the bare interactions and 
have no such a problem. 

We will show in the following that, even after the E-DMFT approximation, the effects 
from these A terms still cancel exactly in all measurable quantities and only the physical 
interactions determine the physics. In doing so, we also establish the equivalence of our 
E-DMFT with the existing ones. The physical reason behind this perfect cancellation is 
that the static quartic interactions, including the A terms we introduced, do not get any 
renormalization in DMFT, even in the presence of those non-local interactions. 

To proceed, we need to check the effects on the phonon dynamical Weiss functions of 
the arbitrary A dependence we added in Eq. @j for ensuring the positive-definiteness. These 
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are easily seen by checking the high frequency behavior of the phonon Weiss functions. In 
this hmit, the phonon self-energy goes to zero at the rate of the inverse frequency square, 
since it is approximately proportional to the local two particle Green's function as shown in 
Eq.(pOD (remember V approaches a finite constant in the same limit). So from Eq.(p!9|) we 
have: 

D'^liiuJn) "-^ _ ^ Vk,o.p = -A„^ (21) 

k 

and thus 

Vapiii^n) -Aa/3 (22) 

We can make a shift of the dynamical phonon Weiss functions by defining (in the matrix 
form) 

V{zun) = V{lUn) + A (23) 

The Weiss function defined above approaches zero in the high frequency limit. We then 
integrate out the auxiliary fields in the effection action Eq.(|13|). After rearranging the 
variables and introducing a new set of the auxiliary field 0o.o5 the effective action becomes: 

Sf = -fdT f dr'Y.ciAr)Q-Ur " r')c,Ar') -\ f dr f dr' ^ 
Jo Jo 2 Jo Jo ^^^^0 

: (poair) : Pj(r - r') : 0o/3(r') : +U dr : no,T(r) :: no,|(r) : - / rfr ^ : 0oa(^) - 5'o,a(r) : 

Jo Jo 

(24) 

It can be shown that, while the electronic Green's functions remain the same, the new local 
phonon Green's function D^°'^ is related to the old one by 

(25) 

One can also find the relation between the phonon self-energies: 
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U-\iuJn) = Il-\iuJn) + X. (26) 

The self-consistent conditions, Eqs.(0) and (|T9|), are transformed to 

f-(z^„) = + Yl-\tuj^)]-\ (27) 

k 

D'^'^itujn) = - E[-^fc"' + n(2cu„)]-\ (28) 

k 

The self-consistent loop can be formed in the same way as before. 

So we have shown that the arbitrariness in the choice of A does not affect the physical 
quantities. This point is best seen through Eq.(^) where a simultaneous shift of Vk and 
cancels out exactly. The quantities related to the auxiliary fields, including those described 
in Eqs. (|^), (H), (126|) , and (pSD , do depend on A. But as we pointed out, they are not 
quantities experimentally measurable. 

The formalism shown above, Eqs.(p^)-(p8D, is equivalent to, although more generic than, 
those discussed in the existing literatures, due to (i) the more general form of the interac- 
tion we have taken, (ii) the consideration paid to possible broken symmetry, and (iii) the 
introduction of the A terms ensuring the positive-definiteness of the interaction matrix. The 
method allows us to study the models with general interactions, like the antiferromagnetic 
spin exchange. 



Before we leave this part, it should be pointed out that the formalism we just devel- 
oped can readily be applied to the general case containing electron-phonon and long range 
Coulomb interactions. There is basically only one change needed. Due to the dynamics of 
the real phonons from themselves, we have an additional term (iun)'^ in those related equa- 
tions, including (|T0|), ([TT|), ([TSl), (|T9|), and (]28|). The electron-phonon identities still hold in 
the new formalism. Without any further change, the continuous auxiliary boson fields can 
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be used to describe the real phonon fields within the same formalism. 

IV. APPLICATION I: PHASE DIAGRAM OF THE 3D U-V MODEL AT 

HALF-FILLING VIA E-DMFT 

As an example of the practical implementation, in the following we apply the E-DMFT 
formalism based on the continuous auxiliary field approach to a 3D Hubbard model with a 
nearest neighbor density- density repulsion (the U-V model). We will be interested in the 
case at half-filling. This is a much simplified version of the model we have investigated 
in the last section. Especially, there is only one local shift, A = Aqo, needed to make the 
interaction matrix positive-definite. Our purpose is to demonstrate the implementation of 
E-DMFT and investigate the behavior of the density response function in approaching the 
CDW phase transition. Under the given condition, this model Hamiltonian allows three 
different phases: the Mott insulating phase when U is dominant; the band insulating phase 
with CDW when V prevails; the metallic Fermi liquid phase when the kinetic energy over- 
comes the interactions. In the MI phase the system can develop an antiferromagnetic long 
order if the magnetic frustration is weak enough. For simplicity, in the following we will 
consider the system at a high enough temperature such that the MI phase is paramagnetic 
due to the strong thermal fluctuations. A finite temperature phase diagram ai j3 = 5.0 is 
presented in Fig.||. In the temperature region where our study is performed, the Mott tran- 
sition is actually a crossoverB so the phase diagram in Fig.|l| should be viewed as a qualitative 
representation of the actual phase diagram at T = 0. 

The 3D U-V model allows us to illustrate the novel aspects of our methodology, including 
the implementation of E-DMFT with an auxiliary field, carrying out QMC simulation with 
repulsive interactions, cancellation of the A dependence, and investigation of the density- 
density response. Physically, we have in mind the following questions: (1) How is the MI-FL 
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transition affected by the non-local interaction V7 (2) How does the charge density insta- 
bility develop in the E-DMFT equations? (3) How does a frequency dependent effective 
interaction affect the quasiparticle properties? These questions can not be addressed in the 
simple DMFT studies where the V interaction is handled at the Hartree level. 



The following is the 3D model Hamiltonian we are going to study: 
H = IK^iiCcja + h.c.) +UY,-nir--nii- 

i,j o" i 

+ VY,{: hi^ : + : nn :)(: hj^ : + : hj^ :), (29) 

where represents a pair of nearest-neighboring sites. In this special case, the Fourier 

transformed off-site interaction is given by 

Vfc = V[cos{k^) + cos{ky) + cos(fc^)]. (30) 

We perform QMC simulation similar to that in Ref.[0]. The fermionic part of the impurity 
model is handled by the standard Hirsch-Fye algorithm^. The statistical weight from the 
part of the continuous bosonic fields is obtained directly by computing the corresponding 
Boltzmann factor. We use here a semi-circular density of state for the electronic degrees 
of freedomlii. The bandwidth is set to be the energy unit. We consider paramagnetic solu- 
tions at a finite temperature. Since the system is exactly at half-filling, we use in QMC the 
particle-hole symmetry accompanied by the reversal of the phonon displacement to increase 
the efficiency of the simulation. In QMC, we take (3 = 5.0 and At = 1/4. Correspondingly 
the inverse temperature range [0, (3] is devided into L = (3/ At = 20 slices. The typical num- 
ber of sweeps for QMC measurement is 10^ by which all the quantities converge within the 
statistical errors. We actually experience critical slowing down in both QMC simulation and 
the DMFT self-consistent iterations as we approach the critical point around U = 3.0 and 
V = 1.6 (see Fig.|l|). In this region the typical number of iterations needed for convergence 
increases from ten to twenty and we use in QMC measurement 2 ~ 4 x 10^ sweeps. 
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Before we present the results, we want to show that at the numerical level, the A depend- 
ing term does not affect the physical results. We have shown in the previous section that 
in the formulation of E-DMFT the A dependence does not show up in physical quantities. 
However it can happen that if the formalism is so sensitive to the dependence such that in 
a practical calculation like QMC there is always only partial cancellation and a significant 
A dependence remains in the physical measurables. It is also possible that a negligible A 
dependence in the results of the impurity model, which is always there unless one can solve 
the problem exactly, may get magnified during the E-DMFT iterations. In the following 
we present the results calculated at f/ = 3.0 and V = 1.6 which show no A dependence. 
From the phase diagram Fig.|l| one can see that this is the point we have reached closest in 
the parameter space to the finite temperature critical point (CP). If the suggested scenarios 
of the A dependence in the physical results may happen, they most likely happen around 
the CP where the system becomes very sensitive to the extra change from A. In Fig.|] we 
plot the imaginary part of the electron Green's functions calculated at four different A's. 
The result shows no A dependence within the accuracy of the calculationill. Similar for the 
electron density Green's functions as shown in Fig.^. So in both the formulation and the 
practical calculation, the A term does not play any physical role and we can use it safely. In 
all the following calculation, we set A = 2.0 in the QMC simulation. 



A. Near the Mott Crossover Line 

We first present the result showing the crossover between the MI and FL phases. In 
Fig.^ we show the data of Im G^°'^{m / (3) v.s. U aX V = 0.0 and 1.0. The reason we choose 
to plot this function is that it is knownS that the asymptotic behavior at the low frequency 
limit of the imaginary part of the electron Green's function reflects the density of states 
(DOS) near the Fermi surface. For the metallic phase the DOS is finite at the Fermi surface 
and hence Im G{ipn) Const. ^ 0. For the insulating phase, on the other hand, the 
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corresponding DOS is zero and Im G{ipn) -—^ 0. So from the plotting of Im G at the first 
fermionic Matsubara frequency we can study the transition or crossover behavior between 
the two phases. Since we solve the problem at a temperature {(3 = 5.0) higher than the 
critical temperature {(3c ~ 1/0.04, see Ref.P]) of the Mott transition, what we see in Fig.^ 
is that for fixed V, — Im G^°^{m/I3) decreases smoothly as U is increased. This behavior 
represents the crossover between the two phase. There is a small but finite shift in between 
the two curves calculated at = 0.0 and 1.0. This is a result of the competition between U 
and V . In the case oiV = 1.0 the system enters the BI phase for U ^ 2.2 which one can see 
approximately from Figj^. This region is out of the scope of the current approach applicable 
only to a homogeneous system. In the range 2.2 U 4.0, which is, roughly speaking, 
the crossover region, the value of Im G^°^{iT[/ {3) at the same U is always bigger for V = 1.0 
than that for V = 0.0. This means that the existence of the finite V makes the system more 
metallic and so the effective U is smaller, as one anticipates from qualitative considerations. 
If U is increased further so that 4.0 ^ U, both the systems enter the paramagnetic insulating 
phase with literally no difference. 

We then investigate the behavior of the various Green's functions as the transition to- 
wards the band insulating phase is approached by changing V at fixed U. We first show 
the results at f/ = 3.0. We solve the problem for seven different values of the interaction: 
V = 0.5, 1.0, 1.3, 1.4, 1.5, 1.6, and 1.7. In all the results presented below the data from 
the last case are not shown because we already encounter the instability at which the con- 
vergence of the self-consistent interaction is lost. So we can view V = 1.7 as an upper 
boundary of the metallic phase for U = 3.0. In Fig.^ we show the imaginary part of the 
electron Green's function and in Fig.|^ the imaginary part of the electron self-energy. It 
can be seen clearly that, while the trends of the change of the plotted quantities are in the 
direction towards a more metallic phase as V is increased (that is, bigger ImG and smaller 
self-energy at the first several Matsubara frequencies), the magnitudes of the changes are 
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very limited, especially in comparison with the phonon Green's functionsE^ shown in Fig.|^. 
The local phonon Green's function, which is related to the density susceptibility as we have 
shown in the last section, increases significantly as we approach the phase transition point. 



As a result of screening the phonon self-energy is negative. Moreover, |n(ici;„)| is a mono- 
tonically decreasing function of the frequency, since the screening becomes less effective 
as the frequency increases. Hence the instability, if it may happen, will first show up at 
the wave vector k = q (tt, tt, tt) and the frequency ci;„ = cuq = 0, where the product 
VqliiioJo) = —3VIl{iuQ) has the biggest positive value. In Fig.H we show the plotting of 
—3VIl{iuJo) v.s. V at the given U. The trend is obvious for the product to approach "1" 
where the corresponding denominator in Eq.(|3lD vanishes. This establishes the picture that 
as the transition point is approached, the denominator disappears first at q and ujq, which 
corresponds to an instability against the homogeneous ground state with a static CDW at 
wave vector (7r,7r,7r). This is a typical phenomenon in the Green's function description of 
phase transitionsii, although the quantities involved here are non-perturbative. The insta- 
bility is signaled by a frequency of phonon becoming neg ative§'0. At this point an ordered 
mean field state would be the correct solution. It should be noted that, even when the 
instability happens, one may still continue the paramagnetic solution of E-DMFT by taking 
the principle part in Eq. (|3TD . But in 3D the convergence near transition is impossible be- 
cause the derivative of the phonon Green's function with respect to the control parameter 
V becomes infinite. 



We can analyze the behavior from the self-consistent equation (0), which can be rewritten 



as 




(31) 
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B. Frequency Dependent On-Site Interaction 



Another interesting property to investigate is the effective on-site density-density inter- 
action, which is defined at A = 

Ues{iuJn) = U + V{iuJn) 

and is frequency dependent. From Fig.^ we can see clearly that as the transition is ap- 
proached, there is a tendency of the softening of the effective interaction at the zero fre- 
quency. In our model the frequency dependence of the local U(,q is due to screening of the 
bare interaction by the intersite Coulomb interaction V. Notice however that a frequency 
dependent U occurs more generally in realistic models of correlated electron due to intrasite 
screening by other local orbitals as a recent local GW calculation0 shows. This is the first 
E-DMFT study of a model where the interaction U is frequency dependent. We should 
stress that this behavior of the single particle Green's function can not be described by an 
ordinary DMFT with fixed U . In Fig.|lO| we show how the frequency depending effective U 
changes the single electron behavior. We plot in Fig.|10| the imaginary part of the electron 
self-energy as a function of the Matsubara frequency. At high frequencies, the self-energy 
from E-DMFT coincides with that calculated using the bare U alone. This tells that the 
screening effect is not effective in the high frequency limit, same as we can see from the effec- 
tive U plotted in Fig.p. In the low frequencies, the self-energy deviates to that of a smaller 
effective U . In Fig.|l^ we plotted the self-energies calculated at = and U"^^ = Ucs{i^^n) 
evaluated at the lowest and the next lowest Matsubara frequencies. We can see that at 
the first Matsubara frequency the E-DMFT self-energy is closest to that given by Ucs{i^^i) 
and V = 0. We can understand the situation by thinking that there are two different t/'s 
controlling low and high frequency regions separately. Some effective and screened U is in 
charge of the low frequency behavior while the bare U works in the high frequency region. 
In between there is a kind of the crossover connecting the two. The results shown here 
suggests that a frequency independent effective U is not enough to capture the physics in 
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the entire frequency range. 

The set of diagrams presented above are plotted very close to the line of the Mott 
crossover (see Fig.Q). In the following, we show two other sets of data which are plotted in 
the metallic and the Mott insulator phases, respectively. 



C. Metallic Phase 



First, we show the results at f/ = 2.0 and increasing V in Figs.|TT]-|l3|. As can be seen 
from Fig.|I|, at this set of parameters the system is in the correlated metallic phase. To be 
concise, we show here three representative plottings, that is the electron and phonon Green's 
functions and the electron self-energy. One can see that in this phase the change of the sin- 
gle electron Green's function is very limited as the transition is approached at V ~ 0.95. 
Meanwhile, the single electron self-energy changes quite a lot in the low frequency region, 
showing the stronger cancellation between the effects from U and V and thus the more 
significant reduction of the self-energy as V is increased. The change of the phonon Green's 
function is again much bigger than those of the electrons. 



D. Paramagnetic Mott Insulating Phase 



The second case with U = 4.0 is shown in Figs.|Tj-|T6[ Here we work in the Mott in- 
sulating phase. One can see from the asymptotic behavior of the Green's functions in the 
low frequency limit that in entire the range, especially near the transition at Vc ~ 3.4, the 
system is still in the Mott insulating phase. Meanwhile, the corresponding phonon Green's 



functions plotted in Fig. 16 change a lot. 
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E. Phase Diagram 



Finally, by literally sweeping across the U-V space, we are able to establish the finite 
temperature phase diagram presented in Fig.|l|. The phase transition from the metallic and 
the MI phases to the BI phase is determined unambiguously from the breaking down of the 
convergence of the E-DMFT iterations. We locate the crossover line between the FL and 
MI phases by search the points of (U,V) at which ImG(zpo) = —0.5. While its specific value 
is arguable, this criterion works well practically in the sense that right around the CP it 
suggests, {Uc, Vc) ~ (3.0, 1.6), we experienced the strongest critical slowing down. 

Two remarks are in place concerning the qualitative features contained in this finite tem- 
perature phase diagram. First, the slopes of the boundaries of the FL phase are positive 
on both sides. This actually refiects the competition between U and V: The existence of 
a finite and small V requires a bigger U in order to make the Mott transition or crossover. 
Similarly a finite and small U makes it harder to develop the CDW. Second, the effects of a 
finite V is much stronger than that of U, because of the coordination number, which is six 
in the current case. The above features will retain in the phase diagram at T = 0. 

V. E-DMFT PLUS C-DMFT FOR SYSTEMS WITH TWO SUBLATTICES 

Next we consider E-DMFT on systems with two interpenetrating sublattices. The cur- 
rent study is useful in the situation when the two sublattices are not equivalent in the sense 
that, while it is homogeneous within each of them, the order parameter is different in the 
two sublattices. One then needs in E-DMFT a cluster containing at least two neighboring 
sites. It is interesting, though, to notice that the formalism we are going to develop also 
applies to the homogeneous systems. In this case the cluster plays the role to improve the 
description of the spatial correlations. The formalism described in this section can be easily 



21 



extended to problems where clusters of bigger sizes are needed. E-DMFT was combined 
with Dynamical Cluster Approximation (DCA) in Ref.P]. 



Under the given conditions, the nearest neighbor hoppings and interactions alway connect 
the two different sublattices, the next nearest neighbor ones are within the same sublattice, 
etc. We rely on the external magnetic fields introduced at the beginning, Eq.(|ll), to lift 
any possible degeneracy in the ground state. To illustrate the basic idea while avoid any 
unnecessary repetition (as it will turn out, the E-DMFT with two sublattices shares many 
properties with that for the homogeneous system) , we work on a model with only nearest 
neighbor hoppings and density- density interactions, besides the on-site energy (the chemical 
potential) and the Hubbard interaction. The Hamiltonian reads: 

H = -t i^\i aCBj,a + h-C.) - ^ fiXi,anXi,a 

{Ai,Bj)a Xi,a 

+ UY^ nxi,^nxi,i + J2 ('^Ai,] + nAi,i)VAi,BjinBj,^ + nBj,i), (32) 

Xi (Ai,Bj) 

where every site is label by Xi with X labeling the two sublattices, X = A,B, and i the 
coordinate within the sublattice. {Ai, Bj) represents a pair of neighboring sites. We choose 
the chemical potential consistent with the two sublattice picture: 

Ua. X = A 

fJ'Xi,a = < (33 j 

[fiBa X = B 

We set for nearest neighbors V(Ai,Bj) = V 0. We then introduce the A term same as before 

Vxi,Yj = ^^XY^ij — Vxi,Yj 

with A a constant which ensures the V2X2 niatrix to be positive-definite. After normal 
ordering the operators in the interactions and perform the Hubbard-Stratonovich transform, 
we have the effective action: 

/■/? 

S = driY, [cxiA^)dTCxiA^) - ^^f,anxiA^)] 
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{Ai,Bj),a 



Xi 



+ ^ E : Mr) : V^ly^ : 0y,-(T) : - : (/.^^(t) :: [nxiA^) + nxiA^)] (34) 



with 



fJ'X,a = - -J2V{Xi,Xj)i[^XjA^) + ^Xj,li^)])^ 

^ j 

= U + A, 

where X = B ii X = A and vice versa. The Green's functions we are going to use are 
defined as follows: 



Gr(zr|zV)'^J:^-(T.cx,.(r)c^,,,(r')), 



t (^'\ 



X^^(.rKV) -(T, : [nx,T(r) +nxa(r)] :: [nw,T(r') + nwa(r')] :), 
The Dyson equations are now 2x2 matrix equations: 



(35) 
(36) 
(37) 





^AB 


-1 

(A;, ipn) = 


iPn + A'Ao- -tk 




"SAA 


SAB 


{k,ipn). 


(38) 








-t-k iPn + I^Ba 




5.BA 









-1 






nAA nAB 




nBA 



{k, iun) (39) 



D^A ^AB ^ _2^^ 

(A;, ia;„) = - 
DBA ^BB _2y_^ X 

In the above equations and for all those with two sublattices, we always define the momentum 
in the reduced Brillouin zone. If the lattice under consideration is of supercubic type in 
d-dimension, one can easily find that for the nearest neighbor hopping and interaction, 
tk = tZ]f=iCos/cj and Vk = V J2i=i cos ki. Similar as that for the homogeneous system, we 
can derive an identity relating the phonon and the electron density Green's functions 
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A. A. 





A 


~2Vk 




nAA 


{k,iujn) = 




+ 






-2V-k 


A 




nBA 



(40) 



We are now in the position to introduce the E-DMFT approximation. Following the 
same procedure as before, we can write down the effective E-DMFT action by using the 
cavity construction. What we do here is that we first integrate out all, including those in 
both sublattices, but two neighboring lattice sites, one from each of the two sublattices. In 
this way, we keep a cluster containing two representative lattice sites. The cluster plays the 
role as a composite impurity which is coupled to the self-consistent fermionic and bosonic 
baths. As we have mentioned earlier, it is foundU that the Hartree terms from the non-local 
interaction across the cluster boundary contribute to the effective action. However, since 
the Hamiltonian we use here is prepared in such a way that there is no longer Hartree terms 
contained in the interaction, the procedure towards E-DMFT becomes very straightforward. 
The effective action is given by: 



3^'' = - f dr f dr' clAr)[Gl 

Jo Jo xY,<x 

1 rP 
-o/ dr rfr'^:0x(r):[D^ 
Z JO Jo ^ 

-t- f drY,: nx,T(r) :: nxa(r) : - f^r^ Mr) :: KtW + Tix.iW] : (41) 
Jo X "^^ X 

with X, Y summed over A, B. From the effective action, we can measure the impurity Green's 
functions and calculate the self-energies by using the local Dyson equations. The self- 
consistency is reached by identifying the impurity Green's functions with the local Green's 
functions which are given as follows: 



^loc,AA 


^locAB 




^AA 


^AB 








^loc,BA 


(^loc,BB 


k 









t-k ipn + /iBa 



(42) 
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^loc,AA £)loc,AB 
£^loc,BA £^loc,BB 



E 



^AA ^AB 
^BA ^BB 



(k, iuJri] 



E 



A -2Vk 

-2V^k A 



+ 



(43) 



^loc,AA ^loc,AB 



^loc,BA ^loc,BB 



E 

k 



^ ,A.A. ^ ,A.B 
X X 



^ ,BA ^ ,BB 
X X 



E 

k 



A 

-2l^-fc A 



+ 



(44) 



Same as before, the A dependences cancel out exactly. Combining the last two of the self- 
consistent equations, (^) and (0), we obtain the following identity: 



X 



locBA, 



[ 


1 


+ 




1 





(45) 



It is not difficult to check that the above equation can be obtained directly from the effec- 
tive action, Eq. (PD . One needs only to write down the phonon Green's function and then 
integrate out the auxiliary fields. By using the local phonon Dyson equation, which is the 
matrix version of Eq. (|16D , one recovers Eq. (|4^) . 



This completes the formulation of E-DMFT with two sublattices. One can see that the 
theory easily combines E-DMFT with C-DMFT. We employ here a cluster of two sites, with 
the application in mind which will be discussed in § |V li| . There is, however, no difficulty to 
extend the formalism to clusters of any size. 
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VI. GW METHOD COMBINED WITH E-DMFT 

As we have mentioned, C- allows to pick out a representative lattice cluster, in- 

stead of a single site, in order to describe a many body system. This makes possible to treat 
the finite range interaction as well as the broken symmetry phase within the cluster. The 
advantage of C-DMFT is that it solves exactly the cluster so that the spatially non-local 
correlations within the cluster are automatically taken into account. In combining with 
E-DMFT, C-DMFT is also able to handle interactions with range beyond the cluster size, 
as we have shown in the last section. However, the price one has to pay is that in solving a 
cluster, a lot more technical resources are needed. 



In this section we propose a less computationally intensive prescription as compared to 
the E-DMFT -|- C-DMFT procedure. It is based on the following physical idea. In real 
materials, the on-site Hubbard interaction U is much larger than the non-local ones. Hence 
the local interaction has to be treated non-perturbatively (namely with DMFT) in order 
to obtain the local self-energy. Meanwhile it is legitimate to make perturbative expansion 
to obtain the non-local part of the self-energy in the spirit of the GW method^. The 
original GW method computes a screened Coulomb line W by summing Random Phase 
Approximation (RPA) diagrams and obtains the one electron self-energy by considering the 
lowest order graph in W, hence the name GW. The E-DMFT-GW approach is derivable 
via the Baym-Kadanoff functional^. The functional derivatives of the 2-particle irreducible 
part of the Baym-Kadanoff functional, ^{G,D), with respect to the full Green's function 
give the corresponding self-energies. ^{G,D) is constructed with the full Green's functions 
G and D and the interaction vertices (As discussed in §VIII| , the choice of the phonon field 
should be done judiciously). The E-DMFT-GW method consists of approximating $ (see 
Fig.p!7D by the leading order non-local graphs and evaluating the rest of the functional $ in 
the local approximation. The E-DMFT-GW self-energies are given by 
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S^G,D) ^E-DMFT 



n..(G, D) = ^ 5.n^r°^"^(G, + (1 - 5.)njf (G, d) (46) 

For the approach to be derivable from a functional, S(G', -D) and n(G', D) has to be calculated 
self-consistently. The GW method has been applied in ah initio calculations of semiconduc- 
tors since the original works by G. Strinati, H.J. Mattausch, and W. However, 
it was pointed out that in LDA-GW for electron systems^ the self-consistency results in 
incorrect one-electron spectra and that it is better to compute S(Go) with Gq the unper- 
turbed Green's function instead. Better total energy, though, is obtained from 
We believe that our proposal resolves this contradiction. The E-DMFT iteration obtains 
the largest self-energy term (the on-site one) self-consistently and non-perturbatively. The 
GW approximation is used for the smaller term (the off-diagonal one). In our model cal- 
culation we find the difference between E-DMFT + non-self-consistent GW and E-DMFT 
-|- self-consistent GW is small. This can be generalized in a straightforward way to realistic 
multiband situations. 



We need to discuss more specifically the non-local self-energy diagrams in our generalized 
GW approach in combination with the E-DMFT. We identify two such contributions. The 
first is the boson exchange diagram which is of the same form as that in the GW method 
[Fig.[l8|(Bl)]. We require that the two vertices Fs be local and come from different lattice 
sites, giving rise to the off-diagonal self-energy. Fig.^(Bl) uses the full electron-phonon 
vertices, instead of the bare ones, in the exchange diagram of the self-energy. Those local 
vertices, which can be measured in E-DMFT-QMC, are defined through the following Green's 
functionEi: 

G'"\T,a,,T^a2-^n) (T,c.,(ri)4(r2)0(r3)) 

= f d<Yl fdr^Y. f dT',G'"^{na,\T[a[)G'"%r^a^\r^a'^) 
Jo , Jo , Jo 
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rJr^(r^a;,r^a^;r^)D'°=(r^|r3). (47) 

Unlike the skeleton diagram commonly used for the electron self-energyil where one of the 
vertices should be bare to avoid over-counting, the diagram Fig.|l8](Bl) uses two full local 
vertices. Our requirement, that the two vertices be from different lattice sites, ensures that 
there is no over-counting. While the skeleton construction uses one bare and one full vertex 
to produce the exact self-energy, our method uses two full local vertices to produce the lead- 
ing non-local correction. We should also remark that both the electron and phonon lines 
appeared in Fig.^(Bl) represent the non-perturbative Green's functions from E-DMFT. 
Especially, the phonon Green's function plays the role as the screened Coulomb interaction 
which in the original GW is obtained by using RPA. 



There is another contribution, which originates from the local interaction U at the second 
order [Fig.|l^(B2)]. The reason this contribution is important is that usually U is much bigger 
than V. The effective local vertex F^"'^ is defined in the following way: 

G(2)'i-(riai,r2a2|r3a3,r4a4) = (T.c.,(ri)c,,(r2)4(r3)4(r4)) 
= Gi°^(riai|r4(T4)G'°^(r2a2|r3a3) - G'''''{T2a2\T,a,)G'°'inai\r,as) 



+ fdr[Y. fdr'.Y. fdr'.Y. f dT',Y.G'^\r,a,\T[a[)G'^%T,a,\T',a^^ 
Jo , Jo , Jo I Jo , 



TY{T[a[, TW2\'r'A. rIa;)G^°^(r^a;|r3a3)G"°^(r:a;|r4a4) (48) 

If there is no external magnetic field, only two spin configurations are allowed in the two 
particle Green's function as well as the vertex, that with all the spins in the same direction 
and that with two spins up and two spins down0. Once we get the local vertex, the corre- 



sponding contribution to the off-site self-energy can be constructed as shown in Fig.l8(B2), 



with again the vertices coming from different sites and ensuring no double counting of the 
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diagrams. It should also be pointed out that the diagrams contained in Fig.0(Bl) and (B2) 
are totally different. This can be seen easily by comparing the non-local lines in the two 
diagrams. 

To compare the importance of the two terms, we can investigate their scaling behavior 
with the spatial dimensionS. We should keep in mind that both G and W scale as 1/ \fd for 
nearest neighbors in real space, with d the dimension, and both the interaction vertices are 
local which means they do not scale. We then see that the leading electron-phonon contri- 
bution scales as and that of the on-site interaction as l/d'^/^. In the infinite dimension 
limit there is no doubt that the electron-phonon contribution is more important. However, 
if we work in finite spatial dimensions (usually < 3) and since the on-site interaction V is 
likely much bigger than the off-site one, V ^ these two can be of the same order practically. 
This actually happens in the example in 1-D which we will show in the next section. 

In the same spirit one can also obtain the leading non-local correction to the phonon 
self-energy, i.e. Fig.|18[B(3), which is the leading non-local correction in terms of the E- 
DMFT interactions. Due to the two electron lines, the diagram scales as l/ii when the two 
contributing lattice sites are nearest neighbors. 

In practical calculations, we first solve E-DMFT iteratively and obtain all the local self- 
energies and the interaction vertices. We then apply GW approximation to calculate the 
non-local self-energies by using the Green's functions obtained from E-DMFT. As in the 
original GW, we assume here that the corrections do not change dramatically the physical 
properties of the system so that we are allowed to use the quantities from E-DMFT directly. 
The last step is to use the approximate self-energies in the exact Dyson equation so that all 
the Green's functions, with both spatial and temporal dependences, can be calculated. 

The GW contributions to the self-energies are calculated in the real space-time as follows 
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for i ^ j (see Fig.0) 



Jo Jo I Jo Jo 

rP rP rfS 1-13 1-13 /-/3 

-/ dnY. dT[Y. dr.Y. dr^Y. dr^Y. dr'.Y 
Jo ^0 ^, Jo Jo ^, Jo ,3 Jo ^, 

(^0-, ri(Ti|r2CT2, T-ia3)G{iTiai\jT[a[)G{iT2a2\jT2a2)G{jT3a'r^\iT3a3)T^l''j{^^^ ^-aaalr^ai, rV) 

(49) 

n^^(.rljy) = ['dnY fdr[Y fdr^Y [' dr^Y 
J^ T:^ Jo Jo Jo ~7 



r!,-(riai, raaa; r)G(^riai|jr{ai)G(jr^a^Kr2a2)r!,-(r^a^, r{a;; r') (50) 

In the above equations we labeled the vertices by the lattice site index with in mind the 
possibility of inequivalent sublattices. There is a symmetry one can use in the calculation 
at i 7^ j: 

G^{it\jt') = [G.(jr'|zr)]t (51) 

It is both physically transparent and technically convenient to perform the generalized GW 
calculation in the coordinate space and imaginary time@. It is also very easy to extend the 
expressions to the systems with different sublattices, as we will show in the next section. 



VII. APPLICATION II: ID BAND INSULATOR VIA E-DMFT PLUS GW 



Since we want to investigate if the GW method can improve the E-DMFT results, we 
need to know the corresponding exact solution of the model under investigation. In this 
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section, we implement the E-DMFT of two sublattices for a ID U-V model with an alter- 
nating chemical potential. This model can be solved exactly at T = via Density Matrix 
Renormalization Group (DMRG)0S. The model is relevant in the study of the interplay be- 
tween the electronic correlation and the electron-phonon coupling in the mixed stack organic 
compoundJ^i'S and the ferroelectric perovskites0. The phase diagram has been studied^. 



Because we solve the impurity model in E-DMFT using QMC simulation which works 
at finite temperatures, we need to make the comparison in the band insulating phase where 
all the excitations are gapped so that the effects of the thermal excitations are suppressed 
at low enough temperatures. We also want to make the comparison in the parameter space 
where the quantum fluctuation is strong enough so that the standard mean field solution 
does not work. It is known that in such a case one needs C-DMFT of at least two sites to 
get good agreement with the exact solution^. Our goal here, though, is mainly to see if the 
GW method can improve the E-DMFT results. 



Due to the above reasons, the E-DMFT we are going to use in this section is slightly 
different from that described in §|V| which combined E-DMFT with C-DMFT. The pure E- 
DMFT (without C-DMFT) for two sublattices is established in the following way. We first 
choose a representative site from sublattice A, integrate out all the other sites in both the 
sublattices, and obtain an effective impurity action for this site. Then, from the nearest 
neighbors of this site, we choose another representative site, which obviously belongs to sub- 
lattice B, and repeat the same procedure. The two effective actions reached in this way are 
the same as that for a homogeneous system we described in ^|ITT[ They are formally inde- 
pendent to each other at the level of the impurity model. Of course the two are connected, 
at the self-consistency, through the Dyson equations which are the same as those given by 
Eqs. ( |38| ) and (|39D, except the off-diagonal self-energies are now zero. Technically, one can 
easily understand the structure of this E-DMFT by imposing the requirement on all the 
corresponding equations in §0 that the impurity model be restrictively local. Then all the 
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off-diagonal dynamical Weiss fields, and thus the off-diagonal self-energies, are zero. How- 
ever, the impurity Green's functions still have non-local contributions, as is evident from 
Eqs.p^)-(P^. This scenario of implementing (E-) DMFT, in midway between the single 
site (E-) DMFT and the cluster one, has been used successfully in treating systems with 
inequivalent sublattices while avoiding the heavy calculations needed in C- DMFT0. 

We study here the following Hamiltonian: 

H = -t^(c|,aCi+i,<7 + h.c.) -^fJ.ini,a 

i,cr i,tj 

1 1 

i i 

We consider the special case with alternating chemical potential: 

fJ^i = (53) 

In this case we know the exact forms of the hopping matrix element and the non-local 
interaction. We can write down for the off-diagonal terms in the Dyson equations (|38|) and 

m 

tk = te'^ cos k, Vk = Ve'^ cos k (54) 

The momentum k is again restricted in the reduced Brillouin zone, given by — 7r/2 < k < 7t/2 
in the current case. The phase factors exp{±ik) come from the fact that each site (A,i) 
[(B,i)] has two neighbors, one within the same unit cell, (B,i) [(A,i)] and the other comes 
from the cell to the left (right), (B,i-1) [(A,i+1)]. Hence the later contributes a momentum 
dependent phase factor. Remember when we use the above two equations for self-consistency 
in E-DMFT, the off-diagonal self-energies for both electrons and phonons should be set to 
zero due the assumption of locality: 

Sf(zp„)=0, H^^(za;„) = 0. (55) 
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Because of the dimensionality, the momentum summations needed in calculating the local 
Green's functions can be carried out exactly, which give: 

with 



and 



k 

4y2nxx 



(57) 



In the above, we used again the notation that X =B if X=A and vice versa. As we have 
noted, the solution of the impurity model in the current case consists of two independent 
parts, one for each representative lattice site, and each of them are exactly the same as that 
for a homogeneous system. The only difference comes in at the self-consistent conditions 
which are given by the above pair of equations. 



After we get the solution of the impurity model, either within a single iteration or after 
the convergence of the E-DMFT iterations, we can perform the GW perturbative calcula- 
tions. To illustrate the idea and see qualitatively how the GW sclf-cncrgy can improve the 
results, we consider here the simplest and the most important GW contributions, those from 
the nearest neighbors. They contribute directly to the off-diagonal self-energies, while for 
the diagonal ones, we use those from the E-DMFT calculation. All the other contributions 
are neglected because the Green's functions decay exponentially as the spatial separation 
increases in the band insulating phase. 



We can now write down the expression for the self-energy matrices as follows: 
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5.AB 


[k, ip„) ^ 

\ " lit/ 


^DMFT.AAj'^ J 


j^GW.AB/- J ifc L 


(58) 
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5^GW,BA(^^p^)g-*fc pQg^ 


SDMFT,BB(,p^) 




nAA 


nAB 




nD^^^T,AA(^^^) 


nGW,AB(^^^)gifc COsfc 














(59) 


nBA 


nBB 




nGW,BA(^-^^)g-ifc pQg^ 


nDMFT,BB(^^^) 





The momentum dependences of the off-diagonal terms in the above two equations come in 
for the same reason as those in the free electron and phonon propagators. In the current 
situation both the GW terms we discussed in the previous section contribute as given by 
Eqs.(|49|) and (|50|). To make life easier, we make a further approximation which replaces the 
full interaction vertices by their bare values: 

r!j°'=(riai, Tsaa; rg) = b^^^^^bijx - T^bijx - rs), (60) 

T^ri'Ti^^i^ T-aCTalrgCTa, na^) = US^^^^J^^^^^S{ti - T2)5{ti - Ts)S{ti - n). (61) 
Then Eqs. (^91) and (0) are greatly simplified and give: 

E^W(^r|jrO = -G.{^T\Jr')Di^T\JT') - U'G^.i^r\Jr')G^,iJr'\^r)G^{^T\JT'), (62) 



U'''^ = Y.GMr\Jr')G.{Jr'\^r). (63) 

a 

In reaching the numerical results we are going to present, we use A = 2.0 for the positive- 
definiteness of the effective interaction matrix. In every iteration, the impurity model is 
solved via QMC by 10^ sweeps. To reach the E-DMFT convergence, ten iterations are usu- 
ally needed. 



We first show in Fig.|T^ the temperature dependence of the imaginary part of the elec- 
tron Green's function calculated at f/ = 5.0, V = 0.5, and /i = 2.0. The three cases are 
calculated at three different inverse temperatures, (3 = 5.0, 8.0, 10.0, with the corresponding 
At = 0.25, 0.20, 0.25, respectively. It is obvious that the three sets of the data lie on a 
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single smooth curve. Actually this same situation happens for all the other quantities we 
measured which are not shown here. All these suggest that, at the given temperatures, with 
the highest at T = 1/5.0, the thermal fluctuations are already suppressed due the band 
gap and we need not worry about the temperature effects. In the following, we present the 
results calculated at /5 = 5.0. 

In Fig.^ we present the data of the imaginary part of the electron Green's functions cal- 
culated using E-DMFT alone, using GW with the electron-phonon vertex [corresponding to 
Fig.p!8[B(l)], and with the local Hubbard vertex [corresponding to Fig.|l8]B(2)], respectively. 
In the latter two the GW calculations are performed after E-DMFT convergence. The ex- 
act result and the Hartree mean field (MF) result are also plotted as references. From the 
results, one can see that the two terms in the GW correction are of the same order. 



In Fig.^ we show the results of the imaginary part of the local electron Green's function 
from GW calculation after E-DMFT convergence and those using GW within the E-DMFT 
iteration loop. One can see that the difference is very small. The corresponding real part 
is plotted in Fig.|2^. In Fig.^ we show the plotting of the Green's function between a pair 
of neighboring sites. We can also compare the result on the average energy per site e. The 
result of GW after E-DMFT gives e = -1.76 and that for GW within E-DMFT e = -1.74. 
The DMRG finds the exact average ground state energy per site to be e = —2.09. The 
difference between the two GW-I-DMFT procedures is again very minor. From the given 
results and those performed at the other parameter points which are not shown here, we 
conclude that the two procedures, one with GW after E-DMFT and the other with GW in 
E-DMFT iteration, give very close results, although it seems the former is a little better. 



The physical information contained in Figs. p0| - p3| can be understood as follows, (i) In 
one dimension, all the allowed modes of the low energy excitations are those bosonic particle 
hole pairs carrying momentum ~ and A;~2A;i?(=7rat half- filling) with respect to the 
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Hartree ground state. This explains why a C-DMFT calculation with a cluster of only two 
sites gives quite good result^ while the E-DMFT we employed here does not work very well 
at low frequencies. The difference is basically that a model of a single site can only capture 
those modes /c ~ 0; but a cluster of two sites is already good enough for those at /c ~ vr. 
(ii) Since the classical Hartree energy gap, which is given by —U/2 + 2V + fi = 0.5 is quite 
small in this case and both the interactions in Eq. (|5^) are (marginally) relevant with respect 
to the metallic Gaussian fixed point, as we go to lower energy scales and thus longer wave 
lengths, the energy gap gets renormalized significantly. This explains why the exact DMRG 
result is so different from the MF result at low frequencies. The high frequency behavior, 
on the other hand, can be captured fairly well even by the Hartree approximation, which is 
evident from Figs. |20| - p3| . This is the region where all the different approaches converge to 
give the same result, (iii) From the figures of the local electron Green's function one can see 
that at not too low frequencies (basically, those beyond the first two Matsubara frequencies), 
the &DMFT result is much better than the MF result and closer to the exact ones. This is 
consistent with the scaling picture since the contributions to the results at those frequencies 
higher than the gap energy, which is of the order of "1" in the current case, can only be from 
the local behavior and are described fairly well by the E-DMFT. On the other hand, this 
same reason explains the big deviation of the E-DMFT result at the first two Matsubara 
frequencies: They are affected more strongly by those quantum/thermal fluctuations with 
longer wave length, which are mostly neglected by the E-DMFT approximation, (iv) The 
GW works in the way we have anticipated. It incorporates more spatial correlations into 
the self-energy so that the low frequency behavior benefits a lot from the correction. As we 



can see from Figs.^ and ^ the GW contributes to both the real and imaginary parts of 
the local Green's functions corrections of more than 15%. (v) The GW method used here 
has little effect on the nearest neighbor Green's functions because it is designed to improve 
the local Green's functions. 

What we have shown above is that the GW method can be used to improve the E- 
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DMFT results. Due to the dimensionality, the single-site E-DMFT results deviates from 
the exact ones at the lowest frequencies. By incorporating the spatial correlations, the GW 
perturbation contributes a desired correction. Of course in this very case in ID it is known 
that the spatial correlation is so important that a leading order perturbation is not enough 
to recover the exact results. What is important is that the above example shows the GW 
method works in the way as we anticipated. Our major objective is to apply GW-DMFT 
to strongly correlated electronic systems in higher dimensions, we know E-DMFT works 
much better as is evident from the scaling behavior with respect to the dimension. We also 
know that the leading order perturbation in terms of the interaction vertex, the GW, works 
better. We thus have a method which is much easier to handle technically than C-DMFT 
and is able to achieve the same goal to certain extent. 



VIII. FURTHER DEVELOPMENT AND OUTLOOK 

E-DMFT allows to describe the interactions in a more flexible way than we have presented 
so far. The need for such a freedom is evident in the realistic calculation of materials, where, 
instead of the nearest neighbor repulsion considered in this paper we have to treat the 
Coulomb interaction and its multipole expansion. For such a system, the dielectric function 
is given by, in the linear response theory!!: 

e"^ (g, iun) = 1 + Vqxiq, iuJn) (64) 

where Vq is the Coulomb interaction given by, in 3D 

Vq = ^ (65) 

and xilji^n) is the density-density Green's function defined through Eq.(P). One can make 
use of the electron-phonon identity, Eq.(|TO|) (set A = 0) and get: 

e(g, iUn) = 1 - Vgll{q, iUn) (66) 
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Here the phonon self-energy 11 can be understood as the collection of all the electron polar- 
ization diagrams. If one proceeds with the formalism we presented in the previous sections, 
the phonon self-energy 11 within E-DMFT approximation is assumed to be momentum in- 
dependent. However, Eq.(^) is not the correct functional form for an insulator in which 
the polarization should be given by H(g, ~ (ffil^^n) with / weakly q dependent. 



To handle this situation in E-DMFT we have to tailor the formalism to be compatible 
with the functional form of the response function. This leads us to the following generaliza- 
tion of the action discussed in Sections |I| and IITT : 



S — Sn 



Jo „ 



a,b 



(67) 



where we have defined a generalized electron density: 

k Q/3 



(68) 



with a and (3 the spin labels. The indices a and h are used to label the local degrees of freedom 
other than the spin, like the components of the multipole moments. 5*0 is the free action plus 
the local interaction. The key part here is the electron-phonon vertex Aa-a,i3{k + q/2, k — q/2) 
which is momentum dependent. A wise choice of this vertex allows to preserve the physical 
momentum dependence in the response function after the E-DMFT approximation. 

We can define the following Green's function in the matrix form: 



[Diq,T)U = -{TrMq,r)U~Q,0)) 



[x{Q,r)]ab = -{^TPa{q,r)pb{-q,0)) 



Y^Ka,ap{k + g/2, k - q/2)\^.p.{k' - g/2, k' + g/2)x„/3,a'/3'(^, 9, k', -q, 0) (69) 



k,k' 



with 
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Xal3,a'f3'{k, q, T] k' , q' , 0) = -(T^4^q/2,^(r)Cfc_g/2,/3(r)4,^g,/2,„,(0)Cfc'_g72,;3'(0)) 

Same as before, we can derive an electron-phonon identity: 

[Xiq,tc0n)r' = -Do{q) + U-\q,tcUn), (70) 

So far the results are exact. The E-DMFT approximation amounts to mapping the general 
model ( p^ ) to an impurity problem by integrating out all but one lattices site (we consider 
the homogeneous phase here). The only new feature is that the general electron-phonon in- 
teraction vertex is not necessarily local as we had before. The resulting phonon self-energy 
is a function of frequency only. However, the general density-density Green's function x, a 
physical quantity, now contains a non-trivial momentum dependence which is evident from 
Eqs.([69|) and (|70D . Since the electron-phonon vertex A can alway be adjusted by redefining 
the auxiliary phonons, we are able to obtain the desired momentum dependence from the 
physical considerations. 



To see how it works, we go back to the example of the insulator. We need the following 
form of the electron-phonon coupling in order to describe the dipole-dipole interaction: 

Aa-Mk + g/2, k - q/2) = 6^p^ (71) 
with a = x,y, z. Meanwhile the free phonon propagator is of the form 

Do,M = (72) 

in order that the Coulomb interaction be recovered when the auxiliary phonons are inte- 
grated out. Under this interaction vertex, the auxiliary phonon represents an electric field 
mediating the dipole-dipole interaction. Since only the longitudinal field is coupled to the 
dipole moment, we can keep the corresponding phonon mode and discard the transversal 
ones. For the E-DMFT approximation, it is desired to work in the coordinate space where 
the cavity construction is possible. To this end one needs to convert the electron-phonon 
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coupling to real space which can be done by using the Wannier functions. 

After solving the E-DMFT problem, the electron density Green's function is given by 

2 -n(ia;„) , , 

X{q, tuJn) = q 7— j — 5777^ — r (73) 

and the dielectric function becomes 

e{q, iujn) = 1 + 4:7re'^Ii{iuJn), (74) 
which has the correct form for an insulator. 

To conclude, we have introduced in this section a way to tailor the E-DMFT formahsm 
so that the desired momentum dependence can be preserved from physical considerations. 

IX. CONCLUSION 

In the paper we suggested a simple procedure of deriving E-DMFT formalism, that is 
first separating out the Hartree contributions and then making the E-DMFT approximation 
with regard to the fluctuations around the Hartree ground state. This procedure is essential 
in the phase with broken symmetry. It also helps to formulate the C-DMFT. 

We developed an E-DMFT formulation by using a real Hubbard-Stratonovich trans- 
formation. We introduced a local shift to the general non-local interaction to ensure the 
positive-deflniteness of the effective interaction matrix. Our investigation showed that in 
all the physical quantities the effects from the arbitrary shift canceled out exactly. We 
also proved the equivalence of forming the E-DMFT self- consistency by using the auxiliary 
phonon Green's function and the two-electron Green's function. Based on these ideas, we 
derived an E-DMFT of a single impurity site for a homogeneous system with generic two 
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particle interactions. We also presented a formalism of E-DMFT combined with C-DMFT 
for a cluster of two lattice sites, which is generalizable to clusters containing any number of 
sites. 

We suggested a generahzed GW approach to incorporate the spatial correlations into the 
E-DMFT approximation. While the on-site self-energies are obtained non-perturbatively 
through E-DMFT, those relatively weaker off-site contributions can be calculated in a per- 
turbative way. We identified the most important contributions to the non-local self-energies. 

We showed how E-DMFT could be tailored to handle the response functions with non- 
trivial momentum dependence in an insulator. Through the example of the dielectric func- 
tion we exhibited that an appropriately defined electron-phonon vertex was able to keep the 
correct functional form of the response function. 

We implemented a QMC algorithm with shifts in E-DMFT to handle the non-positive- 
definite interactions. This algorithm can be used for a large variety of problems, including 
the Anderson lattice with antiferromagnetic interactions. Two examples of the implemen- 
tations were presented. 

The first example was the apphcation of the single site E-DMFT to the 3D U-V model. 
We studied the behavior of the electron Green's function and the response function as the 
density instability was approached. We studied the crossover between the metallic and the 
Mott insulating phases. We investigated the frequency dependence of the effective on-site 
interaction and showed its impact on the single electron behavior. We showed a finite tem- 
perature phase diagram of the 3D U-V model. 

In the second example, we applied a single site E-DMFT combined with GW method 
to a ID U-V model with an alternating chemical potential. It was found that the GW 
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approach improved the E-DMFT results at low frequencies in the desired direction. We also 
found that in the case under investigation, it made little difference whether or not the GW 
perturbation was performed within the E-DMFT iteration. 

The success of the E-DMFT implementation opens the door to tackle many complicated 
physical problems which could not be handled by simple DMFT or other methods. The 
combination with GW and/or C-DMFT points out a systematic way to improve the (E-) 
DMFT method. 
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FIG. 1. The phase diagram of the 3D U-V model at /? = 5.0. The centers of the symbols 
represent the numerical results whose accuracy is of the order of 0.01. The line bounded the BI 
phase represents the phase transitions from the FL and MI phases. The phase transition line is 
found by approaching the instability of the E-DMFT iteration from the smaller values of V . The 
line in between the FL and MI phases represents a crossover. It is determined by the values of U 
and V at which Im Ga{ipo) = —0.5. 
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FIG. 2. The imaginary part of the electron Green's function at four different values of A calcu- 
lated at [/ = 3.0 and V = 1.6 as a function of the Matsubara frequency0. 
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FIG. 3. The density susceptibility at four different values of A calculated at U = 3.0 and 
V = 1.6 as a function of the imaginary time. 
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FIG. 4. The imaginary part of the electron Green's function at the first Matsubara frequency, 
Im G^°^ {iTT / P) , as a function of the on-site interaction U for two different values of V. 
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FIG. 5. The imaginary part of the local electron Matsubara Green's function for different 
values of the interaction V at fixed U = 3.0. 
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FIG. 6. The frequency dependence of the imaginary part of the electron self-energy So-(ipn) 
at different values of the interaction V at fixed U = 3.0. 
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FIG. 7. The local phonon Matsubara Green's function for different values of the interaction V 
at fixed U = S.M. 
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FIG. 8. The rescaled phonon self-energy at the most instable point with k = (vr, 7r,7r) and 
iOn = for U = 3.0. The diamonds represent the calculated data. The curve is plotted using a 
polynomial fitting of the phonon self-energy 11 (which is even in V) up to V^. From the extrap- 
olation the transition is at ~ 1.85. This value is less accurate than the one from QMC due to 
the non-analyticity near the transition. 
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FIG. 9. The frequency dependence of the effective on-site interaction UesiicOn) at different 
values of the interaction V at fixed U = 3.0. 




FIG. 10. The imaginary part of the electron self-energy calculated using E-DMFT at U = 3.0 

and V = 1.5 and at the corresponding = Ucs{ii^n) for n = 0,1, oo with V = 0. We used 
UesiiLOo) = 2.540 and Uesiicoi) = 2.760 from our E-DMFT calculation. Obviously C/eff(iWoo) = 3.0 
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FIG. 11. The imaginary part of the local electron Matsubara Green's function for different 
values of the interaction V at fixed U = 2.0. 
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FIG. 12. The imaginary part of the local electron self-energy for different values of the inter- 
action V at fixed U = 2.0. 
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FIG. 13. The local phonon Matsubara Green's function for different values of the interaction 
V at fixed U = 2.0. 
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FIG. 14. The imaginary part of the local electron Matsubara Green's function for different 
values of the interaction V at fixed U = 4.0. 
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FIG. 15. The imaginary part of the local electron self-energy for different values of the inter- 
action V at fixed U = 4.0. 




FIG. 16. The local phonon Matsubara Green's function for different values of the interaction 
V at fixed U = 4.0. 



E-DMFT-GW ^ J. 

FIG. 17. The 2-particle irreducible functional ^[G,D] in the E-DMFT-GW approach. 
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A. Interaction Vertices 




(1) (2) 
B. Generalized GW Self-Energies at i ^ j 




(1) (2) (3) 

FIG. 18. A. The bare interaction vertices contained in Eqs.(^ and (p^), that is, the local elec- 
tron-phonon interaction A(l) and the local Hubbard interaction A(2). B. The non-local self-energy 
contributions described in Eq.(^9|) [B(l) and B(2)] for electrons and Eq.(|50|) [B(3)] for phonons. 
The strengths of the leading contributions (from the nearest neighbors) for the three diagrams are 
0[V/d], 0[U^/d?/% and 0[l/d], respectively. Since in each of the diagrams we require the vertices 
be from different lattice sites, there is no double counting. 
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FIG. 19. The imaginary part of the local electron Green's function G^^'^'^^{ipn) at U = 5.0, 
V = 0.5, and fi = 2.0. The results at three different inverse temperatures are showiJU. Within the 
accuracy of the calculation, the three sets of data lie on a smooth curve which means the thermal 
effects on the result has already been suppressed due to the band gap. 




FIG. 20. The imaginary part of the local electron Green's function G "'^' {ipn) at U = 5.0, 
V = 0.5, fj, = 2.0, and /? = 5.0. The data labeled as GW(1) comes from the contribution de- 
scribed by FigJl|B(l) and GW(2) from Fig.|l|B(2). It can be seen that GW(1) and GW(2) makes 
corrections at the same order to the E-DMFT result. 
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FIG. 21. The imaginary part of the local electron Green's function G °''^' {ipn) at J7 = 5.0, 
V = 0.5, n = 2.0, and /3 = 5.0. 
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FIG. 22. The real part of the local electron Green's function G^^'''^'^ {ipn) at U = 5.0, V = 0.5, 
II = 2.0, and /? = 5.0. 
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